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For suspensions of permeable particles, the short-time translational and rotational self-diffusion coefficients, 
and collective diffusion and sedimentation coefficients are evaluated theoretically. An individual particle is 
modeled as a uniformly permeable sphere of a given permeability, with the internal solvent flow described by 
the Debye-Bueche-Brinkman equation. The particles are assumed to interact non-hydrodynamically by their 
excluded volumes. The virial expansion of the transport properties in powers of the volume fraction is per- 
formed up to the two-particle level. The first-order virial coefficients corresponding to two-body hydrodynamic 
interactions are evaluated with very high accuracy by the series expansion in inverse powers of the inter-particle 
distance. Results are obtained and discussed for a wide range of the ratio, x, of the particle radius to the hydro- 
dynamic screening length inside a permeable sphere. It is shown that for x > 10, the virial coefficients of the 
transport properties are well-approximated by the hydrodynamic radius (annulus) model developed by us earlier 
for the effective viscosity of porous-particle suspensions. 

PACS numbers: 82.70.Dd, 66.10.cg, 67.10.Jn 
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I. INTRODUCTION 

One of the theoretical methods to analyze transport prop- 
erties in suspensions of interacting colloidal particles is the 
virial expansion in terms of the particle volume fraction cf>. 
For suspensions of non-permeable hard spheres with stick hy- 
drodynamic boundary conditions, virial expansion results for 
short-time properties are known to high numerical precision 
up to the three-particle level, i.e., to quadratic order in <f) f° r 
diffusion and sedimentation coefficients |0 0], and to third 
order in 4> f° r the effective viscosity |0]. The concentration 
range of applicability of these hard-sphere virial expansion 
results in comparison to simulation data has been discussed 
in |01 ■ Our knowledge on virial expansion coefficients of col- 
loidal transport properties is less developed when suspensions 
of solvent-permeable particles are considered. The theoretical 
description of their dynamics is more complicated since one 
needs to account for the solvent flow also inside the particles. 
Permeable particle systems are frequently encountered in soft 
matter science. Prominent examples of practical relevance 
which are the subject of ongoing research, are dendrimers 
10, 1(1, microgel particles 10-131, a large variety of core-shell 
particles with a dry core and an outer porous layer B10l4l4ll . 
fractal aggregates IU5I1 . and star-like polymers of lower func- 
tionality lfl6ll . 

In a series of recent articles lfl7l - l20ll . we have explored the 
generic effect of solvent permeability on the short-time trans- 
port using the model of uniformly permeable colloidal spheres 
with excluded volume interactions. This simple model is 
specified by two parameters only, namely the particle volume 
fraction <fi = (47r/3)na 3 , where n is the number concentra- 
tion and a is the particle radius, and the ratio x of the particle 
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radius to the hydrodynamic penetration depth inside a perme- 
able sphere. Large (low) values of x correspond to weakly 
(strongly) permeable particles. While the model of uniformly 
permeable hard spheres ignores a specific intra-particle struc- 
ture, it is generic in the sense that the hydrodynamic structure 
of more complex porous particles can be approximately ac- 
counted for in terms of a mean permeability. In our related 
previous publications, using a hydrodynamic multipole sim- 
ulation method of a very high accuracy [21], encoded in the 
HYDROMULTIPOLE program package 10], we have calculated 
the short-time translational diffusion properties [ 17lll8ll . and 
the high-frequency viscosity lfl9l,l20ll of the permeable spheres 
model, as functions of <fi and x. These results cover the full 
range of permeabilities, with volume fractions extending up 
to the liquid-solid transition. 

While the simulation results are important for the general 
understanding of permeability effects in concentrated sys- 
tems, for practical use in experimental data evaluation and as 
input in long-time theories, virial expansion results based on 
a rigorous theoretical calculation are still strongly on demand. 
In fact, the knowledge of the leading-order virial coefficients 
can be a good starting point in deriving approximate expres- 
sions for transport properties, which may be applicable at con- 
centrations much higher than those where the original (trun- 
cated) virial expansion result is useful. An example in case is 
provided in our recent derivation ll20ll of a generalized Saito 
expression for the effective high-frequency viscosity 77^ of 
permeable spheres, based on the second-order concentration 
expansion result, i.e. a Huggins coefficient calculation, of this 
property. 

In 11911 . we have performed virial expansion calculations for 
rjoo. We have investigated therein a simplifying hydrodynamic 
radius model (HRM), where a uniformly permeable sphere 
of radius a is described by a spherical annulus particle with 
an inner hydrodynamic radius a e s(x) < a, and unchanged 
excluded-volume radius a. In this annulus model (HRM), the 
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Huggins coefficient describing two-body viscosity contribu- 
tions has been evaluated and shown to be in a remarkably 
good agreement with the precise numerical data for porous 
particles, characterized by a wide range of permeabilities re- 
alized in experimental systems. 

In the present article, the aforementioned theoretical work 
on the virial expansion of the high-frequency viscosity is gen- 
eralized to short-time diffusion properties. In Sec. [II] the 
virial expansion is performed for the translational and rota- 
tional self-diffusion coefficients D t and D r , respectively, the 
sedimentation coefficient K, and the associated collective dif- 
fusion coefficient D c = KD /S(Q). Here, 5*(0) is the small- 
wavenumber limit of the static structure factor, and D is the 
single-particle translational diffusion coefficient. In Sec. HIH 
we provide highly accurate numerical values for the first-order 
(i.e., two-particle) virial coefficients, Xt(x), X r (x), Xc(x) and 
Xk (x), of D t , D r , Dc and K, respectively, in the full range 
of permeabilities. In Sec. lIVI we also recalculate these virial 
coefficients approximately on the basis of the simplifying an- 
nulus model (HRM). In Sec.|V] we conclude that in the range 
x > 10 typical of many permeable particle systems, the virial 
coefficients are well approximated by the annulus model. 

II. THEORY 

We consider a suspension made of a fluid with shear vis- 
cosity 7/0 and identical porous particles of radius a. The fluid 
flow is characterized by Reynolds number Re<< 1. Outside 
the particles, the fluid velocity v and pressure p satisfy the 
Stokes equations ll22ll23ll . 

770 V 2 v(r) - Vp(r) = 

V-v(r) = 0, (1) 

and inside the particles, the Debye-Biiche-Brinkman (DBB) 
equations i24ll25ll . 

770 V 2 t>(r) — i] k 2 [v(r) — Ui(r)] — Vp(r) = 

V-v(r) = 0, (2) 

where is the hydrodynamic penetration depth. The skele- 
ton of the particle i, centered at r,, moves rigidly with the 
local velocity Uj(r) = Uj + u>i x (r — r^), determined by the 
translational and rotational velocities and u>i of the par- 
ticle i, respectively. The fluid velocity and stress tensor are 
continuous across a particle surface. The effect of the parti- 
cle porosity is therefore described by the ratio x of the parti- 
cle radius a to the hydrodynamic screening length k~ x of the 
porous material inside the particle, i.e. 

x = na. (3) 

Owing to linearity of the Stokes and DBB equations and the 
boundary conditions, the particle velocity U, depends linearly 
on the external forces Fj exerted on a particle j. In particular, 
for two interacting spherical particles, 1,2, in the absence of 
external torques and flows, 

Ui = Mii(ri,r 2 )-F 1 + At * t 2 (r 1 ,r 2 )-F 2 (4) 
U 2 = j u| t 1 (ri,r 2 )-F 1 + i uf 2 (r 1 ,r 2 )-r 2 . (5) 



In this paper, the two-particle translational-translational mo- 
bility matrices /Lt*j(ri,r 2 ) are evaluated using the multipole 
method of solving the Stokes an DBB equations ll2lll . The 
cluster expansion of the above mobility matrices reads, 

^(r 1 ,r 2 ) = /i*«5 ij -l+ / x*f ) (r), (6) 
where r = r 2 — ri and 



is the single porous-particle translational mobility. Here Aiq 
is a single porous-particle scattering coefficient B26I1 . given ex- 
plicitly in AppendixlAl For a non-permeable hard sphere with 
stick boundary conditions, A\ s = 3a/2. 

The single particle scattering coefficients A\ a , with I = 
1,2,3,4,... and a = 1,2,3 [26], are essential to perform 
the multipole expansion. They determine the corresponding 
multipoles of the hydrodynamic force density on a particle 
immersed in an ambient flow; examples are Eqs. (|7]i and d24l) . 
The same coefficients Ai a specify also the corresponding mul- 
tipoles of the fluid velocity, reflected (scattered) by a particle 
immersed in a given ambient flow. This is why Ai a are called 
"scattering coefficients". In the multipole approach, differ- 
ences in the internal structure of particles (e.g. solid, liquid, 
gas, porous, core-shell, stick-slip), are fully accounted by dif- 
ferent scattering coefficients. The other parts of the multipole 
algorithm need not to be changed. The scattering coefficients 
are the matrix elements of two single-particle friction oper- 
ators, Zo and Zo, which determine the hydrodynamic force 
density exerted by a given ambient flow on a motionless and 
a freely moving particle, respectively. 

In the multipole expansion method, the two-particle mobil- 
ity /x( 2 )(l, 2) (e.g. translational one, as inEq. ||6), or rotational 
one, as in Eq. d23l ) is expressed in terms of the single-particle 
friction operators, Zo(i) and Zo(i), with i = 1,2, and the 
Green operator G(l, 2). The later relates the flow outgoing 
from particle 2 and incoming on particle 1 . We can write fi^ 
as an infinite scattering series, 

^ - /x Zo(l + GZ )- 1 GZoMo 

= // ZoGZo/x — /x ZoGZoGZo// + .... (8) 

Since the multipole matrix elements of the Green tensor G 
scale as inverse powers of the interparticle distance r, Eq. ([8]) 
corresponds to a power series in 1/r. Truncating the expan- 
sion at order 1/r 1000 , we obtain a very high precision of the 
mobility calculation, actually much higher than needed for 
any practical applications. 

In the present work, we investigate the short-time dynam- 
ics, at time scales t << a 2 / Dq, where 

D* = k B T& (9) 

is the single-particle translational diffusion coefficient, with 
the Boltzmann constant ks and temperature T. On the short- 
time scale, the system is described by the equilibrium particle 
distribution [27]. In the further analysis, we will need only 
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the small-concentration limit go(r) of the pair distribution 
function, where r is the interparticle distance. For particle- 
particle direct interactions described by a pair potential V(r), 
this pair distribution is go(r) = exp(— V{r)/kBT). For non- 
overlapping spheres of radius a, 



9o{r) 



for r < 2a, 

1 for r > 2a. 



(10) 



The first-order terms in the virial expansion of the short-time 
transport coefficients are obtained by averaging the corre- 
sponding two-particle mobility elements. As the result, the 
virial coefficients are obtained as integrals, which involve the 
mobility elements and go (r) . 

The first-order virial expansion of the short-time transla- 
tional self-diffusion coefficient has the form, 



Dt — Dq(1 + X t (f> + 0{4> 2 )). 
The coefficient At is given by the relation ll28l 



A f = 8 



with R = r/2a and 



i(R)J t (R)R 2 dR, 



J t {R) = ^Tv^ 2 \v). 
Mo 



(11) 



(12) 



(13) 



Here, Tr denotes the trace operation. For the sedimentation 
coefficient, one obtains, 



K=l + \ K (t> + 



where 



and 



5a 3 
+8 

J K {R) = \lr 
Mo 
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+ 00 



[go(R) 

g (R)J K (R)R 2 dR, 



(14) 

l]RdR (15) 
(16) 



tt(2) 

Mil 



(r) 



yi2 (2) (r)-T (r) 



In the above expression, 

T (r) 



8irrjor ' 



(17) 



(18) 



is the Oseen tensor and r = r/r. 

The scattering coefficient A12 for a porous particle llibll is 
given explicitly in Appendix lAl For a non-permeable hard 
sphere with the stick boundary conditions, A\% — 5a 3 /2. 

The collective diffusion coefficient is given by 



D C = D'oK/SQ), 



(19) 



where 5(0) is the zero-wavenumber limit of the static struc- 
ture factor, 5(0) = lim 9 ^o S(q). The first-order virial expan- 
sion of Eq. ( [T9i > has the form, 



D c - + Ao0 + •••)• 



(20) 



For the non-overlapping spheres J27|], 

5(0) = 1 -8(/) + 0((t> 2 ). 

In this case, 



Af7 = A 



K 



(21) 



(22) 



The relation <|22) follows from Eqs. (0, (O and j20l . 

We proceed by analyzing the short-time rotational self- 
diffusion coefficient. In the absence of external forces and 
flows, the two-particle rotational-rotational mobility matrices 
/x[J(ri , Y2) satisfy the relation analogical to Eqs. (|4j-(|5j, with 
the translational velocities replaced by the rotational ones, and 
the forces replaced by the torques. The two-particle cluster 
expansion now reads, 



M«"(ri,r 2 ) = HoSijl 



M? 2) (r), 



with 



Mo 



Snr] Ai 



(23) 



(24) 



The scattering coefficient An for a porous particle 1 26] is 
given in Appendix lAl For a non-permeable hard sphere with 
the stick boundary conditions, A^ = a 3 . 

The virial expansion of the rotational self-diffusion coeffi- 
cient is 



where 



and 



with 



D r =D r {l + \ r <j> +...), 



D r = k B Tfig 



(25) 



(26) 



/+00 
g (R)J r (R)R 2 dR, (27) 



J r (R) = ±TrnYi (2) (r). 
Mo 



(28) 



III. RESULTS 

To evaluate the first order virial coefficients A = At , Xk , A r , 
we calculate two-particle mobility matrix elements, perform- 
ing a series expansion in powers of 1/r up to the order 1000, 
as described in Sec.[n] Integration with respect to the particle 
positions has been performed analytically term by term, using 
the expressions given in the previous section. The virial coef- 
ficients A have been evaluated for a wide range of x. Selected 
results are listed in Table|U All displayed digits are significant. 
(More data are available on request.) 

In Fig. [T] the first order viral expansion of the short-time 
translational self-diffusion coefficient of a porous-particle sus- 
pension is compared with the accurate simulation results per- 
formed in Ref. 11711 for volume fractions <f> < 0.45. The first 
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TABLE I: First virial coefficients At, \k and A r for the short- 
time translational self-diffusion, sedimentation and rotational self- 
diffusion, respectively. 



X 


A t 


\k 


A r 


3 


-0.2497 


-3.4451 


-0.03257 


4 


-0.4159 


-4.1066 


-0.06336 


5 


-0.5692 


-4.5539 


-0.09682 


6 


-0.7021 


-4.8722 


-0.12956 


7 


-0.8149 


-5.1084 


-0.16012 


8 


-0.9102 


-5.2898 


-0.18802 


9 


-0.9909 


-5.4328 


-0.21327 


10 


-1.0598 


-5.5480 


-0.23606 


1 1 


-1.1190 


-5.6426 


-0.25662 


13 


-1.2151 


-5.7884 


-0.29208 


16 


-1.3202 


-5.9380 


-0.33426 


1 o 


-1.3 1 JU 


-O.UUbO 


ft ^AQ'"> 


20 


-1.4161 


-6.0662 


-0.37628 


30 


-1.5499 


-6.2335 


-0.44202 


40 


-1.6190 


-6.3149 


-0.48007 


50 


-1.6610 


-6.3628 


-0.50497 


65 


-1.7001 


-6.4064 


-0.52959 


100 


-1.7460 


-6.4563 


-0.56075 


DC 


-1.8315 


-6.5464 


-0.63054 



spheres. For uniformly porous particles, the boundary 
collocation method was applied by Chen and Cai IBfJl 
to evaluate \ K = -3.46, -5.50, -6.23, -6.44 for 
x 2 = 10", with a = 1,2,3,4, respectively. Com- 
paring their results with our very accurate values, A^ = 
-3.5723, -5.5480, -6.2504, -6.4563, we conclude that 
the uncertainty of their results is decreasing from 3% at a = 1 
to 0.3% at a = 4. The accuracy of the boundary collocation 
method is worse at smaller values of x, i.e. for larger perme- 
abilities. 

In Figs. OU the first-order virial coefficients A are plot- 
ted versus the porosity parameter 1/x, for 1/x < 0.1, i.e. 
for x > 10. For 1/x = 0, i.e. for x = oo, the limit 
of a non-permeable hard sphere with radius a is recovered, 
A hs = A(oo). As shown in Fig. [2] for a low permeability, the 

-4.5, 1 1 1 1 1 
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x=30 
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X=oo 
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0.25 



0.35 



0.45 



-6.5- 



0.02 



0.04 



1/x 



0.06 



0.08 



0.1 



FIG. 2: Two-particle sedimentation virial coefficient \k(x). Our 
precise results for porous particles (solid line) are well-approximated 
by the annulus model (dashed line). 

coefficient Xk ~ A 1 ^ + 10/ x is approximately a linear 
function of 1/x. The coefficients At and X r as functions of 
1/x are shown in Figs. [3] and [4] 



FIG. 1: Translational self-diffusion coefficient D t for a suspension 
of porous particles. Symbols connected by splines (solid lines): ac- 
curate simulation results from Ref. [ 17]. Dashed straight lines: first- 
order virial expansion calculated in this work. 



order virial expansion, see Eq. (fTTl i. can be used as an accurate 
approximation of D t in a very wide range of volume fractions, 
even for relatively large values of x (i.e. low permeabilities). 

In contrast, values of the sedimentation coefficient K 
differ significantly from the first-order virial estimation al- 
ready at rather small volume fractions, see Ref. d. 
Our values of A^(oo) reproduce with a higher accuracy 
the classic Batchelor's result B29I1 for non-permeable hard 
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FIG. 3: Two-particle translational self-diffusion virial coefficient 
Xt(x). Our precise results for porous particles (solid line) are well- 
approximated by the annulus model (dashed line). 




0.02 



0.04 0.06 

1/x 



0.08 



0.1 



FIG. 4: Two-particle rotational self-diffusion virial coefficient 
X r (x). Our precise results for porous particles (solid line) are well- 
approximated by the annulus model (dashed line). 



IV. ANNULUS MODEL 



In the annulus model H31H . a particle suspended in a viscous 
fluid is characterized by two radii, a< and a>. Its hydrody- 
namic interactions are governed by the smaller radius a< . In 
addition, there exist also direct pair interactions. Two parti- 
cles cannot come too close to each other, with the no-overlap 
radius equal to a>. For such a model, the first-order virial 
expansion of transport coefficients has been performed with 
respect to the volume fraction = (47r/3)net^,. The corre- 



sponding first-order virial coefficients X A 
been evaluated as functions of e, where 



- \ A> \ A 



, X A have 



(29) 



The method used to determine X A is described in AppendixlBl 
The calculated values are listed in Table ITR 



TABLE II: First-order virial coefficients Xf, A# and for the 
short-time translational self-diffusion, sedimentation and rotational 
self-diffusion, respectively, for a suspension of the annulus particles. 



c 


A t 




A r 


0.00 


-1.8315 


-6.5464 


-0.63055 


0.01 


-1.7523 


-6.4601 


-0.56666 


0.02 


-1.6793 


-6.3769 


-0.51671 


0.03 


-1.6109 


-6.2962 


-0.47417 


0.04 


-1.5466 


-6.2179 


-0.43699 


0.05 


-1.4860 


-6.1419 


-0.40402 


0.06 


-1.4286 


-6.0680 


-0.37451 


0.07 


-1.3743 


-5.9962 


-0.34791 


0.08 


-1.3228 


-5.9263 


-0.32381 


0.09 


-1.2739 


-5.8582 


-0.30189 


0.10 


-1.2274 


-5.7918 


-0.28187 


0.11 


-1.1832 


-5.7272 


-0.26354 


0.13 


-1.1008 


-5.6027 


-0.23122 


0.18 


-0.9253 


-5.3166 


-0.16974 


0.24 


-0.7595 


-5.0135 


-0.12034 


0.31 


-0.6111 


-4.7051 


-0.08296 


0.45 


-0.4093 


-4.1986 


-0.04242 


0.66 


-0.2401 


-3.6278 


-0.01775 



Now we are going to compare the first-order virial coeffi- 
cients, calculated in the previous section for porous particles, 
with the corresponding results, obtained in this section for the 
annulus model, called also the hydrodynamic radius model 
(HRM). A similar comparison has been done in Ref. B19I1 for 
the effective viscosity. The key concept in this procedure is 
the hydrodynamic radius of a porous particle. For the transla- 
tional diffusion (self-diffusion and sedimentation), the hydro- 
dynamic radius a* ff is obtained from the single-particle trans- 
lational diffusion coefficient, with the use of the relation, 



D 1 



k B T 



(30) 



The dependence of a* ff on the porosity parameter x follows 
from Eqs. (|7]l and ([gj, which determine the translational dif- 
fusion coefficient of a single porous particle ll24ll25l[32ll . and 
Eq. (IA1I ), which specifies the scattering coefficient Aiq. Ex- 
plicitly, 



Oeff 



(,:) 



2x 2 {x — tanh(a;)) 



(31) 



2a; 3 + 3(>-tanri(:r))' 

For the rotational self-diffusion, the hydrodynamic radius 
cig ff is obtained from the rotational diffusion coefficient of a 
single porous particle i32Ll33ll . with the use of the relation, 

k B T 



87"7o« ff ) 3 



(32) 



6 



The dependence of o£ ff on the porosity parameter x follows 
from Eqs. d24l) and (l26T l. with the scattering coefficient An 
given by Eq. (IA2b . As the result, 



'eft 



(,') 



1 + 



3 coth(a;) 



1/3 



For large x, 



'ell 



(x) = at s (x) + O (1/x 2 ) . 



(33) 



(34) 



A porous particle of radius a and the porosity parameter x 
is modeled as an annulus particle, see Fig. [5] by matching its 
geometrical radius a to the annulus no-overlap radius a> = a. 
The smaller annulus radius a< = a e ff is determined by the 
effective hydrodynamic radius a e ft, given in Eqs. d3~TTl-(l3"3"l>. 













a \ 
















POROUS 



ANNULUS 



FIG. 5: The annulus (or hydrodynamic radius) model of a porous 
particle. 

In this way, the annulus parameter e, defined in Eq. d29l ), 
becomes the following function of x, 



e(x) 



a eS (x) 



a ef f(x) 



(35) 



with a e ff(a;) determined by Eqs. (IBTt and (l33l l. 

In Figs. |2]|4j the annulus coefficients A j4 (e), with e(x) spec- 
ified by Eq. d35l l (dashed lines) are compared to the porous- 
particle virial coefficients X(x) (solid lines). For the sedimen- 
tation coefficient, the annulus model is accurate, with a half- 
percent relative accuracy already at x = 20 and a reasonable 
3% precision at x = 5. For the translational self-diffusion, the 
annulus model is less accurate, but still it gives only a 2% er- 
ror for x — 20, and a 5% error for x = 10, and a 7% error for 
x = 5. The least accurate is the annulus prediction for the ro- 
tational self-diffusion, with a 5% error for x — 20 and a 1 1% 
error for a; = 10. Summarizing, the annulus model (HRM) ap- 
proximates well the first virial coefficients of porous particles 
suspensions, in the range of intermediate and small particle 
permeability (i.e. for moderate and large values of x). 



V. CONCLUSIONS 

In this paper, the short-time diffusion properties of dilute 
suspensions of uniformly porous spherical particles have been 
investigated. The first-order virial coefficients A of the dif- 
fusion and sedimentation coefficients have been evaluated as 
functions of permeability. Values of A are well-approximated 
by the annulus (hydrodynamic radius) model, if the parame- 
ter x is sufficiently large, i.e. the permeability is sufficiently 



low. Systematically, the annulus approximation slightly un- 
derestimates the virial coefficients of porous particles suspen- 
sions. For rotational diffusion, a reasonable 5% accuracy of 
this model is reached at x > 20. For translational diffusion 
(collective and self), the comparable or even better 3-5% pre- 
cision is obtained already for x > 10. For the sedimenta- 
tion coefficient, the accuracy is even higher (a 3% precision 
already at x = 5), owing to much larger absolute values of 
A /< = 8 + Ac- 

The annulus model is expected to work well also at larger 
volume fractions, if the porosity parameter x is sufficiently 
large. The accuracy of this approximation at larger volume 
fractions will be investigated in a separate publication. More- 
over, the annulus model is also useful to estimate diffusion 
and shear viscosity coefficients of suspensions made of parti- 
cles characterized by a different internal structure, such as the 
core-shell. This problem will be the subject of a future study. 



Appendix A: Scattering coefficients 

The scattering coefficients for a uniformly permeable 
sphere with x = na, where a is the sphere radius and 1/n 
is the hydrodynamic penetration depth, have the form 13211 . 



A m = 

An = 
An = 

Bl2 = 



(2l+l)gi(x) 



m+i(x) a 2i+i 



l(2l-l)(2l + l)gi(x) 
(l + l)x 2 gi- 2 (x) 



,21-1 



2Z + 3 2(2m)(2Z + 3) 



(Al) 
(A2) 



21-1 {l + l)x 2 

2(21 - 1)(2/+1)' 



a 2 A 
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21+5 
21-1 



,21 + 1 



(I + l)x 2 



a Ai 2 - a 



21+3 



(A3) 
(A4) 



where I = 1, 2, ... and gi(x) = ^Jtt/2xIi +1 /2{x) is the modi- 
fied spherical Bessel function of the first kind. 



Appendix B: The annulus (hydrodynamic radius) model 

For a suspension of particles described by the annulus (hy- 
drodynamic radius) model, the virial coefficients \ A are func- 
tions of the model parameter e, defined by Eq. d29| i and listed 
in Table [TT] In this Appendix, we explain how these values 
have been evaluated. 

The first-order virial expansion of D t / Dq, Dc / Dq, K and 
D r /DQ, can be performed with the use of </>< = (47r/3)na^, 
or 0> defined by the analogical expression, 

1 + X A (e)cb > + O (4>l) - l + \ A (e)<j )< +0{<jP < ). (Bl) 



Therefore, 



A A (e) 



(l + e)s 



(B2) 
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To evaluate A A (e), we now introduce the dimensionless in- 
terparticle distance as R = r/2a < , and we replace the pair 
distribution function from Eq. ( TTOb by the following expres- 



sion, 



9o(R) 



for R < 1 + e, 

1 for R > 1 + e. 



(B3) 



which corresponds to the no-overlap condition at a larger ra- 
dius a>. Then, we apply the Eqs. ( TT2l >. (TToT l and (|27| |, taken in 
the non-permeable hard-sphere limit, .t = oo. We obtain the 
following expressions, 



(e) = 8 



+ 00 



J t (R)R 2 dR 



J t {R)R 2 dR, 



(B4) 



A&(e) 



4^"^ (1+e)2 

/•+oo 

8 / J K (R)R 2 dR 

Jl+e 



= A* - 8 



l+e 



J K {R)R 2 dR, 



(B5) 



/■+oo 

A^(e) = 8 / J r (R)R 2 dR 

Jl+e 



\ hs 



l+e 



J r {R)R 2 dR, 



(B6) 



where 



Mo 



Tr 



Mfi (2) (R) 



^f)(R)- 



(B7) 



In this Appendix, all the mobility coefficients, the associated 
functions J, and the superscript hs , refer to the hard spheres 
with the stick boundary conditions. 

The above formulas have been applied to compute the func- 
tions A/(e), Xi(e) and A;5(e), listed in Table ED The func- 
tions A^(e), Aj$(e) and A^(e) have been also evaluated and 



listed in Ref. I3411 . where they have been applied to model dif- 
fusion and rheology of particles with a hard solid core and a 
thin porous shell. 
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